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We investigate the AC Josephson current through a quantum dot with strong Coulomb interac- 
tion attached to two superconducting and one normal lead. To this end, we perform a perturbation 
expansion in the tunneling couplings within a diagrammatic real-time technique. The AC Josephson 
current is connected to the reduced density matrix elements that describe superconducting correla- 
tions induced on the quantum dot via proximity effect. We analyze the dependence of the A C signal 
on the level position of the quantum dot, the charging energy, and the applied bias voltages. 
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I. INTRODUCTION 

Josephson junctions can be formed by linking two su- 
perconductors via an insulator, a normal conductor or 
a constriction in an otherwise continuous superconduct- 
ing material* 1 ^ Advancements in nanofabrication en- 
abled to contact superconductors with quantum dots 
(QDs), which can be formed in carbon nanotubesj^r— 
in InAs nanowires^— in graphene^ or by means 
of self-organization in InAs with Al electrodes*^— 
One motivation to investigate hybrid superconductor-QD 
device a 17 i 18 is the possibility to tune their properties via 
external electrodes that shift the discrete energy levels 
of the QD. Another feature characteristic for QDs is the 
charging energy that may give rise to effects based on the 
interplay of Coulomb repulsion and superconducting cor- 
relations. It is, therefore, an interesting question to ask 
how the discrete level spectrum and the charging energy 
affect the DC and AC Josephson transport between two 
superconductors coupled via a QD. 

In the absence of a bias voltage, a finite DC current 
can be sustained in such a S-QD-S by the DC Joseph- 
son effect. This has been confirmed experimentally;^ 
which shows that two electrons forming a Cooper pair 
can tunnel coherently one by one through a strongly 
interacting quantum dot. When neglecting the charg- 
ing energy, the DC Josephson effect in such a S-QD-S 
system can be studied within a scattering approach^ 
But also Coulomb-interaction effects have been in- 
cluded in various formalisms as perturbation expan- 
sions in the tunneling HamiltonianSSr— and in the 
Coulomb repulsion ; 24 ' 25 a mean-field approach^ quan- 
tum Monte-Carlo simulations^ a renormalization-group 
technique^ or numerical diagonalization of an effective 
dot Hamiltoniani^Sr— 

A finite bias voltage gives rise to a more complicated 
transport behavior. In addition to a finite DC current, 
sustained by quasiparticle tunneling and (multiple) An- 
dreev reflection, there is a time-dependent component 
due to the AC Josephson effect. Theoretical works for 
this regime have mainly concentrated on the limit of 
vanishing or weak Coulomb repulsion. The DC com- 



ponent has been studied by focusing on single quasi- 
particle tunneling^ 2 - - — using a slave boson mean-field 
approximation^^ or performing a perturbation expan- 
sion in the charging energy^ Multiple Andrecv reflec- 
tion processes not only give rise to a stationary current 
but also lead to higher harmonics contributing to the AC 
Josephson transport. A quantitative description includ- 
ing this interplay has been investigated in quantum point 
contacts by means of a scatterin g 37 ' 38 or a Hamiltonian 
approach^ In a noninteracting quantum dot the depen- 
dence of the different harmonics on the bias voltage has 
been studied^ Also dephasing effects introduced by a 
third, normal electrode added to the S-QD-S setup have 
been investigated within a Keldysh formalism applied to 
a noninteracting system to find that for gradually in- 
creasing coupling to the normal conductor the A C signal 
decreases 4i Further works deal with the time evolution 
of the current after switching on a finite bias voltage ; 42 ' 43 
or with polaronic effects due to coupling to vibrational 
modest 



The aim of this paper is to analyze the AC Joseph- 
son effect through QDs with strong Coulomb repulsion 
that cannot be neglected or treated perturbatively. To 
this end, we extend a real-time diagrammatic approach 
for the DC current presented in Refs. I22U23I to A C trans- 
port. We apply this formalism to a three-terminal ge- 
ometry consisting of a strongly interacting quantum dot, 
which is weakly tunnel coupled to two superconductors 
and one normal conductor, see Fig.Q] and perform a per- 
turbation expansion in tunnel couplings to lowest order. 
To sustain (DC and AG) Josephson currents through the 
QD to first order in the tunnel-coupling strengths, su- 
perconducting correlations must be induced on the QD 
by the proximity of the superconducting leads. We find 
that a large AC current between QD and one supercon- 
ducting lead requires a large proximization of the QD by 
the other superconducting lead which can be achieved by 
tuning the gate voltage accordingly. We discuss the am- 
plitude of the A C components of the current between the 
superconductors as a function of gate and bias voltage. 
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FIG. 1: (Color online) A sketch of the device under inves- 
tigation: a quantum dot with strong Coulomb interaction is 
weakly tunnel coupled to one normal and two superconduct- 
ing leads. 



II. MODEL 

The considered system consists of a quantum dot tun- 
nel coupled to one normal and two superconducting 
leads. Its total Hamiltonian is given by H = i?dot + 
Htun + Y^Hr- The index r refers to the leads and can 

r 

take the values 51, 52, N. The quantum dot is assumed 
to accommodate one spin-degenetate level e. It is de- 
scribed by the Anderson impurity model, 



(1) 



,(t) 



where n a = d\d a is the number operator and d a " an- 
nihilates (creates) an electron with spin a on the dot. 
Coulomb interacting is accounted for by the charging en- 
ergy U for double occupation. The Hilbcrt space of the 
isolated dot is four dimensional and is spanned by the 
kets {\x}} with \ = 0> t> -Jo corresponding respectively 
to empty, singly occupied with spin up, singly occupied 
with spin down, and doubly occupied dot. 

The two superconducting leads are modeled by the 
mean-field BCS Hamiltonian 



H 



BCS 



k,a \ fc 



Cr — k^Crk^ 
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Here, c rk<J is the creation operator of an electron in lead 
r with momentum k and spin a, and A r is the super- 
conducting pair potential. The operator S r ^ annihilates 
(creates) a Cooper pair in lead r, which ensures particle 
conservation. By making use of a Bogoliubov transfor- 
mation, the Hamiltonian can be diagonalizcd, 



E rkll ka lrk<j + ^ r N r 

rka 



(3) 



with the quasiparticle operators 7 r fe CT and the correspond- 
ing eigenenergies E r k = \J (ife — yU r ) 2 + |A r | 2 . The chem- 
ical potential of superconductor r is given by fj, r and N r 



is the total number of electrons, which is the number 
of quasiparticles plus twice the number of Cooper pairs. 
Also the normal conductor can be described by Eq. ([2]) 
by simply setting Ajv = 0. 

Tunneling between dot and leads is described by the 
tunneling Hamiltonian 



-fftnn — 



rka 



t r C. 



H.c. 



(4) 



The tunneling amplitude t r as well as the density of 
states p r are assumed to be energy independent in the 
window relevant for transport. Furthermore, we define 
T r = 2ir\t r \ 2 p r . Finally, we set h = 1 and we reinstate it 
in the units used for the figures. 

Due to tunneling, the superconducting leads can 
induce superconducting correlations on the quantum 
doti 22 ' 23 The resonance condition for this proximity effect 
due to superconductor r is that the energy for the doubly 
occupied dot, 2e + U, equals the energy for an empty dot 
plus an extra Cooper pair in the condensate of the super- 
conductor, 2p r . It is, therefore, convenient to introduce 
the detunings S r = 2e + U — 2p r with r = 51, 52. 



III. METHOD 

A real-time diagrammatic approach to DC transport 
through quantum dots tunnel coupled to normal and su- 
perconducting leads, has been introduced in Ref. H^. In 
the following, we briefly review this formalism and extend 
it to describe AC Josephson transport. The system under 
consideration can be divided into three subsystems, the 
dot, the fermionic states of the leads, and the Cooper 
pair condensates in the superconductors. Since we are 
not interested in the fermionic dynamics of the leads, we 
can trace out their degrees of freedom and arrive at a 
reduced density matrix for the remaining part, i.e., the 
dot's degrees of freedom and the Cooper pair conden- 
sates. Its elements are given by P^ 1 = (£i \p Te d (£2) 1 where 
|£) = \xA n si, n S2}) includes the dot state \ = 0,t, I, d 
as well as the number of Cooper pairs in the two super- 
conductors, 77,51 and ns2, measured relative to an arbi- 
trary but fixed reference. The diagonal elements P^ = P| 
give the probability to be in state £. The off-diagonal ele- 
ments Pj^ 1 with describe coherent superpositions. 

The states £1 and £2 in P^ 1 provide more information 
than is needed to study the electric transport. In fact, 
only the differences of the Cooper pair numbers of the 
states and ^2) are important. Moreover, particle 
number conservation sets the constraint that the total 
number of electrons in has to be the same as in \^)- 
For convenience, we define 

msi,ms2 

(5) 

From the definition it follows that the symmetry re- 
lation P£({n S i,n S2 }) = [P£({-nsi,-ns 2 })]* holds. 
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As a consequence of particle conservation, 1232 is a 
unique function of Xii X2, and nsi- It is, therefore, 
enough to keep track of the Cooper-pair number n of 
one lead only. We choose here lead 51, i.e., n = risi, 
and introduce the definitions P x ( n ) = Py ({ n -> ~ n \) an d 
Pq(ti) = Pq({ti, — n — 1}). To inherit the symme- 
try relation P^(n) — [Pq(— n)] , we consistently define 
P°(n) = P°({n,-n + l}). 

Finally, we collect all the nonvanishing elements of 
the reduced density matrix in the vector 7r(n) = 
(P (n) , P t (n) , Pi (n) , P d (n) , P d (n) , P° (n) ) T . Its dynam- 
ics is governed by the generalized master equation 

— 7r(n)(f) + iE n 7r(n)(t) 
dt 

= £ / d^W(n,n')(t,t>(n')(*'). (6) 

" ' -00 

where l the matrix elements of the kernel 

W XlX j (n, n')(t, t') are the transition rates from an 

X2X2 y ' /v ' y 

initial state at time t' described by P*, 1 (n')(t') to a final 
state at time t described by P^ 1 (n) (t) . For the kernel we 
have introduced a notation analogous to the one adopted 
for the reduced density matrix: W***, 1 (nsi,n' si )(t,t') = 

W*^({n sl7 n S2 } 7 {n' sl ,n' s . 2 })(t 7 t'), where the excess 
number of Cooper pairs in the superconductor S2 is 
fixed by particle conservation. We get 7151 + 7132 to be 
equal to for xi — X2, equal to —1 for xi = d, X2 = 0, 
and equal to +1 for xi = 0, xi — d. In a similar way, 
n' S2 is determined in terms of n' si , x'i; an d X?,- The 
only nonvanishing matrix elements of the matrix E ra are 
E n 11 = 2neV, E n g = S S2 + 2neV = S S1 + 2(n + l)eV, 
and £„™ = — P_„qq, where is the voltage drop 
between S2 and SI and it reads — eV = fis2 — (J-si with 
e > 0. 

The tunneling current between the dot and lead r is 
given by 

t 

I r (t)=eJ2 I dt'e T W r (0,n')(t,t')7T(n')(t'), (7) 

-00 

where e T = (1,1,1,1,0,0) and n = in W r (0, n')(t, f) 
ensure that the final state on the right hand side is diag- 
onal both in the dot state and the Cooper pair numbers. 
The current rates W r (n,n')(t,t') are similar to the gen- 
eral rates W(n, n')(t, t') but take into account the elec- 
trons transferred from lead r to the dot. In addition to 
tunneling currents, there are, in general, displacement 
currents due to the formation of image charges when the 
dot occupation varies in time4£ The displacement cur- 
rents do not play any role for the DC part. But even 
for the AC part they drop out for the symmetrized cur- 
rent Is(t) = [Isi(t) — Is2{t)}/2 when choosing the ca- 
pacitances of the tunnel contacts between the dot and 



the two superconductors symmetrically. Therefore, we 
ignore the displacement currents in the following. 

The frequency of the A C Josephson signal is given by 
the energy difference of a Cooper pair being in supercon- 
ductor SI or 52, i.e., by 2((j,si — = 2eV. Therefore, 
we perform a Fourier expansion by making use of A(t) = 

XT=-oo A n e 2ineVt and A n = (1/7) dtA{t)e.- 2ineVt 
with T = 2ir/ (2eV). Within the diagrammatic approach 
the factor exp(— 2ineVi) appearing in the n-th Fourier 
component of the current simply adds a term 2neV to 
the energy difference of the states on the upper and lower 
Keldysh contour. This term can easily be incorporated 
into the energy difference arising from different Cooper 
pair numbers by shifting W r (0, n')(t, t 1 ) and Tv(n')(t') in 
Eq. to W r (n, n' + n)(t, t') and n(n' + n)(t'), respec- 
tively, i.e., only the 0-th Fourier components of W r and 
7r are needed. Performing the remaining time integral, 
we get for the n-th Fourier component of the current 

^ = e]Te T W>y)7r(n'), (8) 

n' 

with W = f_ d£'W(£, t')e l0+t being the zero frequency 
Laplace transformed rate, that does not depend on the 
final time t. The 0-th Fourier components 7r(n') are read- 
ily obtained from the 0-th Fourier component of Eq. ^ , 

iE„7r(n) = W K n')Tv{n') (9) 

n 1 

together with the normalization condition e T n(n) = 
5 n _o- In summary, the n-th Fourier component of the 
current can be evaluated within the diagrammatic tech- 
nique in exactly the same way as the DC current (see 
Ref. HH) but allowing for off-diagonal final Cooper pair 
states (n ^ 0) in W r (n,rt') on the right hand side of 
Eq. ©. The diagrammatic rules to calculate the kernels 
W and W r are given in Appendix [21 

IV. RESULTS 

In the following, we perform a systematic perturba- 
tion expansion of 7r(n), W(n,n'), W r (n,n') and /" in 
the tunnel-coupling strengths, T = maxjrsi, T$2-, Tat}. 
Since we assume the tunnel couplings to be weak, we 
restrict ourselves to lowest (first) order for the kernels 
W(n, n') and W (n, n'). In addition, we concentrate on 
the limit of an infinite superconducting gap in the leads, 
A — > 00, i.e., quasi-particle tunneling between dot and 
the superconductors is suppressed. As a consequence, the 
current into the superconductors is exclusively sustained 
by Cooper pairs. The normal lead affects the occupation 
of the quantum dot, which, in turn, affects Cooper-pair 
transport. Even a weakly tunnel-coupled normal conduc- 
tor influences AC Josephson transport between the two 
superconductors . 

For A -> 00 and to first order in T, all matrix el- 
ements of W(n,n') and W r (n,n') that require either 
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higher-order tunneling or a finite superconducting gap A 
in the leads vanish. The only non-vanishing ones enter- 
ing Eq. ([8]) are readily evaluated, see Appendix [B] This 
results in 



TS 1 =ieT sl [P°(n+1)-P d (n-1)] 



(10) 



i.e., the n-th component to the current into supercon- 
ductor 51 is fully determined by the density matrix ele- 
ments P$(n - lj and P°(n + 1) = [P${-n - 1)]*. The 
latter describe superconducting correlations induced on 
the quantum dot due to the proximity effect. 

The Cooper pair degree of freedom n in 7r(n) intro- 
duces an apparently infinitely large number of density 
matrix elements that are all coupled to each other via 
Eq. (jnj). However, in the limit of a large bias voltage 
as compared to the tunnel-coupling strength, \eV\ T, 
only very few of them need to be taken into account. This 
is a consequence of E„ appearing on the left hand side 
of Eq. (0). Most of its matrix elements are of order eV, 
while W(n, n') on the right hand side scales with V. This 
mismatch defines a hierarchy in powers of T/(eV) for the 
density matrix elements. The lowest order contains all 
matrix elements of 7r(n) for which the corresponding E„ 
is zero or of the order of T. This includes all diagonal 
matrix elements P* (0) (and excludes all elements P* (n) 
with n 7^ 0). The next order contains all matrix elements 
of 7r(n), that can be connected to lowest order ones by 
the kernel W(n, n'). The only off-diagonal matrix ele- 
ments that can be reached from the diagonal ones for 
A ->• oo and to first order in T are Pq(— 1) = [P°(l)]* 

and Pq (0) = [P°(0)]*. If the gate voltage is tuned such 
that the quantum dot is in resonance either with super- 
conductor 51 or 52, namely \5si\ <T or \8s2\ ^ I\ then 
Pq(— 1) or Pq (0) already belong to the lowest order in 
the hierarchy, indicating strong proximity effect with su- 
perconductor 51 or 52, respectively. But in any case, all 
off-diagonal matrix elements except Pq (— 1) = [P° (1)] 
and Pq (0) = [Pj'(O)] can be dropped for describing the 
current into the superconductors. 

Starting with the kinetic equations with off-diagonal 
final states, 

i(6 S2 + 2neV)P d (n) =W$(n, 0)P (0) + W^(n, 0)P d (0) 
+ W d d (n,n)P d (n), (11) 

for n — —1,0, and using the rates listed in Appendix IBl 
we find that the required off-diagonal density matrix el- 
ements are related to the diagonal ones via 



P d (0) = 



2A 



[Po(0) - P rf (0)] (12a) 
[Po(0) - P„(0)] , (12b) 



where we defined the complex resolvents Asi = Ssi + 
iW d d (-l, -1) and A S2 = S S2 + iW d o d (0,0). The ex- 
pressions for Wqq(— 1,— 1) and Wqq(0, 0) are given in 
Appendix [Bj Their imaginary parts can be interpreted 



as the renormalization of the detuning 6si and 832, re- 
spectively, due to the tunnel coupling to the normal 
lead. Their real parts provide a width to the reso- 
nances. For a systematic perturbation expansion, we 
may replace the full expressions for the resolvents A r 
by their leading-order term only. To do so, we need to 
distinguish the two cases of the quantum dot to be on 
or off resonance with superconductor r. On resonance, 
I'M iS r, we find that A r starts to first order, i.e., we 
can omit the S r appearing in the argument of the Fermi 
and the digamma functions to arrive at Re A r = S r + a 



with a" = 
and ImA 



£5 Re 



4' 



2nk B T J V 2 M B T 

T N [1 + f{e) - f(e + U)}. Off resonance, 
|<5 r | 3> r, we can replace the resolvent by its zeroth-order 
term, A r = S r . 

In order to determine the non- vanishing elements of the 
reduced density matrix we also need the kinetic equations 
with diagonal final states, 



°=E w ^( '°) i v(o) 



x' 



2 £ Rete d (0,n')P V) 



(13) 



n'=-l,0 



The rates with diagonal initial and finite states are re- 
lated to single-electron tunneling between dot and normal 
conductor. In contrast, the rates connecting superposi- 
tions between a doubly-occupied and an empty dot to a 
diagonal state require tunneling of one Cooper pair from 
or to the condensate of a superconducting lead. As a 
result, we find 



P (0)-P d (0) = 



l-.f(e)-f(e + U) 



l + f(e)-f(e + U) + Im £ ^ 

r=Sl,S2 

(14) 



Plugging this into Eqs. (|12a[) and (|12b[) and employing 
Eq. (|10l) yields the current into superconductor 51 for 
all values of e. This is what we use to calculate all the 
curves in the figures. The resulting formulae can be sim- 
plified further to obtain compact analytical results after 
specifying whether the quantum dot is in resonance with 
one of the superconductors or not. 

First, we observe from Eq. (fTU)) that only the zeroth 
and first Fourier component of the current flowing into 
superconductor 51 are nonvanishing. For the DC cur- 
rent, Igf = Zgj = 2er s:L ImP d (-l), it is important 
whether the dot is in resonance with the same super- 
conductor the current is measured in. On resonance, 
\Ssi\ < r, we find 



T DC _ 

1 si — 



eT N T 2 sl [1 - /(e) - f(e + U)} 
(S S1 + of + T 2 S1 +T%[1 + /(e) - f(e + U)] 



(15) 

which starts in first order in T. Off resonance, the DC 
current starts only in third order in T, i.e., vanishes to 
the order considered here. In summary, a DC current 
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flows only between the normal lead and the supercon- 
ductor which is in resonance with the quantum dot. We 
explicitly checked that Igf = —1$°, which guarantees 
current conservation. 

We now turn to the AC current Igf(t). It can 
be decomposed into a cos- and a sin-term, Igi(t) = 



0.15 



li?' cos cos(2eVt)+I^ sin sin(2eVi) with I*? cas = I 1 



si 
1st 



2er S iImP o d (0) and Ig[ 



AC.sin 



SI 
(Ik 



SI' 



1st) 



2er S iReP o d (0). The current l£f(t) in 52 can be ob- 
tained from I§i {t) by replacing T^i <h> T S2 and usi ^ 

MS2- 

We immediately see that the behavior of the AC 
tunneling current into one superconductor depends on 
whether the quantum dot is in resonance with the other 
superconducting lead, which supports the interpretation 
that, to lowest order in T, the A C tunneling current into 
superconductor 51 is sustained by oscillations of Cooper 
pairs between lead 51 and the quantum dot that is prox- 
imized by the tunnel coupling to lead 52. 

Similarly as the DC current into lead 51, the cos-part 
of the AC current starts in third order in T as long as the 
quantum dot is off resonance with lead 52. On resonance, 
\o~S2\ ^ r, we obtain 



rAC.COS 

'SI 



eT N T S iT S2 [i-/( e )-/( e + t7)] 

(6s2 + ay + r| 2 + [i + /( £ ) - /( £ + utf ' 

(16) 

The sine-part, on the other hand, is on resonance, \Ss2\ ^ 
r, given by 



AC,sin 
SI 



ersir S2 (SS2 + o-) 



{SS2 

1 - 
X 

1 + 



^S2 



[i + /( £ )-/( e + 



f(e)-f(e + U) 



f(e)-f(e + U) ' 
while off resonance, | <5^2 1 I\ we find 



AC.sin 
SI 



er.siT 



SI 1 S2 



'52 



1-, f(e)- f(e + U) 
1 + f(e)-f(e+U) 



(17) 
(18) 



It is remarkable that the amplitude and phase of the 
AC tunneling current into superconductor 51 does not 
depend on the chemical potential usi but only on the 
detuning Ss2 between quantum dot and the other su- 
perconductor 52. The chemical potential fisi only en- 
ters the oscillation frequency, given by the bias voltage 
2eV = 2(fisi — MS2)- This can be interpreted in the fol- 
lowing way: the tunnel coupling to the superconducting 
leads induces superconducting correlations on the quan- 
tum dot. The total proximity effect is given by the sum 
of the contributions stemming from the two supercon- 
ductors. The AC Joscphson current between quantum 
dot and 51, however, only probes the proximity effect 
induced by superconductor 52, held at a different chem- 
ical potential. 

Another interesting feature of the expression for the 
AC current is the behavior for V — > 0. The ampli- 
tude of the sine first harmonic exactly reproduces the 
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FIG. 2: Amplitude of the AC Josephson current in units of 
2eTs/h as a function of dot level position e/U. The tem- 
perature is UbT = (7/10, the bias voltage to the normal 
lead is hn = 0, and the bias voltages /isi and fis2 ap- 
plied to the superconductors are, respectively, given by a) 
-3*7/2 and 3(7/2, b) -(7/4 and [7/4. The solid curves are for 
Tjv = Ts/10 = (7/100, the dashed lines for T N = T S = (7/10. 
The resonance conditions are indicated by vertical lines. 



DC Josephson transport between the two superconduc- 
tors that was discussed in Ref. This crossover has 
been discussed before in quantum point contact o 38 ' 41 as 
well as in the three-terminal setup under consideration in 
this paper in the noninteracting Hmit4a We also notice 



that the cosine term, / 



ytc.cos 
si 



vanishes for T 



N 



while 



the sine term remains finite. 

In the following, we discuss the amplitude \l£ c \ and 
the phase <fi of only the symmetrized current Ig C '(t) = 
Vsi (*) - J S2 C (*)]/ 2 = \ I s C \ sin(2eW + 0), for which the 
displacement currents drop out. Without loss of general- 
ity, we choose the reference energy such that = 0. 
Furthermore, we concentrate on the limit of symmet- 
ric tunnel coupling to the two superconductors, Tsi = 
Ts2 = Tg. We distinguish the two cases of ~ and 

rv «r s . 

In Figs. [2] and Figs. [3J we plot the amplitudes of the 
AC Josephson currents as a function of the dot level 
position e for five different values of the applied bias 
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voltages usi and Thereby, the voltage is applied 

cither symmetrically or asymmetrically with respect to 
fj> 7v . Coulomb blockade suppresses the AC Josephson ef- 
fect for — U < e < 0, since in this region the quantum 
dot is predominantly singly occupied. This suppression 
is also present for DC Josephson transport through the 
quantum dot, as has been observed experimentally^ and 
discussed theoretically^ Outside this region, an occupa- 
tion of the quantum dot with an even number of electrons 
is possible and, hence, Josephson transport is possible. 

The dot level energies e around which the resonance 
conditions \Ssi\ < r and \Ss2\ T to be in resonance 
with superconductor 51 and 52 are fulfilled are given by 
e « /isi — U/2 and e « fi S2 — U/2, respectively. Around 
these points, which are indicated in the figure by dotted 
lines, the AC tunneling current between quantum dot and 
the other superconductor, 52 and 51, respectively, is en- 
hanced. Depending on the biasing, each of the two reso- 
nances lie either inside or outside the Coulomb-blockade 
gap. In the first case, AC Josephson transport is sup- 
pressed at the resonances and the amplitude of the AC 
Josephson current shows a local maximum close to the 
edge of the Coulomb-blockade gap. In the latter case, a 
maximimum is clearly visible whenever the dot is reso- 
nance with one of the superconductors. 

In Figs. a) and [3] a), both resonances are outside, in 
Fig-Elb) one is outside and one inside, and in Figs. [5]b) 
and [3]c) both are inside the Coulomb-blockade region. 

The phase of the AC Josephson current for symmet- 
rically applied bias voltages is shown in Fig. fJJ We find 
that there is 0-tt transition when crossing through the 
resonances. This means that away from the resonances, 
only the sin-term contributes to the AC Josephson cur- 
rent, while the cos-term vanishes. At the resonance, how- 
ever, the sin-term changes sign, i.e., goes through zero, 
while the cos-term remains finite. For Tn = Ts, dashed 
lines in Figs. [2] and [3] , this 0-7r transition is not visible 
in the amplitude but only in the phase (the phase corre- 
sponding to the asymmetric case, Fig. [31 is not shown). 
In contrast, for <C Ts, solid lines in Figs.[2]and|3j the 
0-7r transition is indicated by a sharp dip in the ampli- 
tude as well. This is a consequence of the fact, that the 
cos-term starts linearly in r^v while the sin-term remains 
finite even for — > 0. 



FIG. 3: Amplitude of the AC Josephson current in units 
of 2eTs/h as a function of dot level position e/U. The tem- 
perature is ksT — (7/10, the bias voltage to the normal lead 
is fix =0, and the bias voltages fisi and ^132 applied to 
the superconductors are, respectively, given by a) —U and 
3(7/2, b) -(7/4 and 3(7/2, and c) -(7/4 and (7/3. The solid 
curves are for Tm = Ts/10 = (7/100, the dashed lines for 
rjv = Ts = (7/10. The resonance conditions are indicated by 
vertical lines. 



V. CONCLUSIONS 

We have analyzed AC Josephson transport through a 
single-level quantum dot tunnel coupled to two super- 
conductors and a normal conductor. The amplitude and 
the phase of the AC Josephson current depend both on 
the gate and bias voltages. As a function of gate volt- 
age, there are two resonances: When the quantum dot 
is in resonance with one superconductor, the AC Joseph- 
son current between dot and the other superconductor is 
enhanced. There is a 0-n transition at each of the res- 
onances. For small tunnel coupling to the normal lead, 
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FIG. 4: Phase of the AC Josephson current in units of n as 
a function of dot level position e/U. The parameters are the 
same as in Fig. [2] a) and b). 



this 0-7T transition is accompanied with a sharp dip in 
the amplitude of the AC Josephson current. Inside the 
Coulomb-blockade region, i.e., the region of gate volt- 
ages for which the quantum dot is predominantly singly 
occupied, the AC Josephson current is exponentially sup- 
pressed. The frequency of the A C Josephson oscillations 
is given by the voltage difference between the two super- 
conductors. 
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Appendix A: Diagrammatic Rules 

for evaluating the generalized rates 



The rules 
W xlxl ({™i' n 2},-!XX 2 }) arc as follows: 



(1) Draw all topologically different diagrams with fixed 
ordering of the vertices in the real axis. The vertices 
are connected in pairs by tunneling lines carrying energy 
LUi. The tunneling lines can be normal or anomalous. 
For each anomalous line choose the direction (forward or 
backward with respect to the Keldysh contour) arbitrar- 



iiy. 

(2) For each vertical cut between two vertices assign a 
factor 1/(AE + irf) with n = + , where AE is the differ- 
ence between the left-going and the right-going energies, 
including the energy of the dot states E x , the tunnel- 
ing lines uii, and the energy difference in Coopcr-pair 
condensates Sep. The latter is increased (decreased) at 
each vertex of an outgoing (incoming) anomalous line at 
which the arrow is opposite to the arbitrarily chosen line 
direction. 

(3) For each tunneling line assign a factor 
±T r D r (cOi)f±(cji), where f+(ui) = f r (ui) = 
[l + exp(w.i -/i r )/(fc B r)] _1 and f~(oJi) = 1- f r {u>i), and 
D r("i) = Tl '"'T'ia i? g (l CJ - Mrl - |A r |). The upper 
(lower) sign applies for lines going backward (forward) 
with respect to the Keldysh contour. For anomalous lines 
multiply an additional factor— ±sign(o;i — fi r ) [J.^l | • 
Moreover, assign a factor e _1 * r for an outgoing and e 1 * 7 " 
for an incoming anomalous line. [For normal leads, only 
normal lines with D r (uii) = 1 appear.] 

(4) Assign an overall prefactor — i. 
Furthermore, assign a factor —1 for each 

a) vertex on the lower propagator; 

b) crossing of tunneling lines; 

c) vertex that connects the doubly-occupied dot state, 
\d) = d|d||0), to spin up, | f); 

d) outgoing (incoming) anomalous tunneling line in 
which the earlier (later) tunnel vertex with respect to 
the Keldysh contour involves a spin-up dot electron. 
[The factors in c) and d) arise due to Fermi statistics from 
the order of the dot and lead operators, respectively] 

(5) For each diagram, integrate over all energies Wj. 
Sum over all diagrams. 



The generalized 
W^' r ({rM,n 2 },K,n' 2 }) 
lowing way: 



are 



current 
evaluated 



rate 



.47 



in the fol- 



(6) Multiply the value of the corresponding general- 
ized rate W**, 1 ({n-i, n%\, {n' l7 n' 2 }) with a factor given by 
adding up the following numbers for each tunneling line 
that is associated with lead r: 

a) for normal lines: 1 if the line is going from the lower 
to the upper, —1 if it is going from the upper to the lower 
propagator, and otherwise; 

b) for anomalous lines: 1 for incoming lines within the 
upper propagator and outgoing lines within the lower 
propagator, —1 for outgoing lines within the upper prop- 
agator and incoming lines within the lower propagator, 
and otherwise. 



Appendix B: Rates and Current Rates 

In this appendix, we list all the rates and current 
rates entering the calculation. Thereby, we omit all rates 
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that can be obtained from the listed ones via the sym- 



metry relations W XlXl (n, n') 
J X2X2 y ' 

W XlX r r (n,n>) - ^ X8 ^' r ' 
X2X2 v ' 



W X2X H-n,-n') 
X1X1 v ' ' 



and 



1. Superconductor 



W xx ,(0,Q) = W xx (0,0) with 



W M {Q, 0)=T N f(e) 
WW(0,0)=IV[l-/(e)] 
W da (0,0) = T N f(e + U) 
W rTd (0,0)=T N [l-f(e + U)], 



To lowest order in T and for A — > 00. many rates and 
currents rates involving a superconducting tunneling line 
vanish. The non- vanishing ones turn out to be indepen- 
dent of the Cooper pair numbers of the condensates. We 
find 

W™(n, n - 1) = W$(n, n + l) = T sl /2 
W d d (n, n-l) = <?(n, n + l) = -T S i/2 
W^(n,n) = W$(n,n) = iT S2 /2 
W${n, n) = W$(n, n) = -ir S2 /2 , 



for the rates and 



IF, 



Od.Sl 



(n,n- 1) = W dd ' S \n,n- 1) = -iT sl /2 



(n, n) = W dd ' S2 (n, n) = -ir S2 /2 , 



for the current rates. 



2. Normal Conductor 

Changing the state of the dot due to tunneling from 
and to the normal conductor is described by the rates 



and Woo(0,0) = -2W CT o(0,0), W CTCT (0,0) = -W 0ct (0,0)- 
W da {0,0), as well as ^(0,0) = -2^(0,0). The 
nonvanishing current rates are given by WzL, (0, 0) = 

W XX ',' N {0,0) with 



wiS(o,o) 

^(0,0) 
W^(0,0) 



-IV[l-/( e )] 
T N f(e + U) 
-T N [l-f(E + U)} 



Finally, we need the kernels W^{n : n) for n = —1,0. 
For n = — 1, we find 



-Re 
-Re 



-T N [i + f( E -S sl )-f(e + U-S S i)] 

\2 2-nk B T 

\2 2irk B T 
For W$(0,0), we get the same but Ssi being replaced 
by S S 2- 



W d d (-l,-l) 

7T 

_ ]£n 
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